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RELATIVISTIC TODA CHAIN AT ROOT OF UNITY II. 
MODIFIED Q-OPERATOR 
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Abstract. Matrix elements of quantum intertwiner as well as the modified Q-operator for 
the quantum relativistic Toda chain at root of unity are constructed explicitly. Modified 
Q-operators make isospectrality transformations of quantum transfer matrices so that the 
classical counterparts of Q-operators correspond to the Backlund isospectrality transforma- 
tions of the classical transfer matrices. Separated vectors for the Functional Bethe Ansatz 
are constructed with the help of modified Q-operators. 
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Introduction 

This paper is the continuation of Ref. |1[] where the quantum relativistic Toda chain at a 
root of unity was investigated. In the previous paper we have established several features of 
the chain, concerning the relationship between its finite dimensional (spin) counterpart and 
its nontrivial classical counterpart. 

For several integrable models, based on the simple Weyl algebra at iV-th root of unity, the 
C-numerical iV-th powers of Weyl generators form a classical discrete integrable system. So 
the parameters of the unitary representation of Weyl algebra form the classical counterpart 
of the finite dimensional integrable system (i.e. spin integrable system), see 0, §, |]. |j. 
Most important finite dimensional operators, arising in the spin systems, have functional 
counterparts, defined as rational mappings in the space of iVth powers of the parameters. So 
the finite dimensional operators are the secondary objects: one has to define first the mapping 
of the parameters, and then the finite dimensional operators are to be constructed in the 
terms of initial and final values of the parameters. Usual finite dimensional integrable models 
correspond to the case when for all the operators involved the initial and final parameters 
coincide. Algebraically these conditions are the origin e.g. of Baxter's curve in the case of the 
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chiral Potts model [H, [7| , or of the spherical triangle parameterization for the Zamolodchikov- 
Bazhanov-Baxter model. 

The simplest case of such combined system is the so-called quantum Toda chain at root 
of unity. It differs significantly from the usual Toda chain || as well as from the relativistic 
Toda chain with arbitrary q |§. For RTC at root of unity L-operators, quantum as well 
as classical, and the quantum intertwiners have rather simple form ||. In [Q it was shown 
that the trace of monodromy of quantum intertwiners, giving usually Baxter's Q-operator, 
has the nontrivial classical counterpart. Under several special additional conditions, the 
classical counterpart of Q operator is the Backlund transformation [|T0| , |TT| for the classical 
counterpart of relativistic Toda chain. Spin counterpart of the corresponding Backlund 
transformation makes an isospectral transformation of the spin system. Parameters of Weyl 
algebra are described in the terms of r-functionF], and the transformation mentioned maps 
(n — l)-solitonic state into n-solitonic state. Initial homogeneous state is supposed to be 
0-solitonic, while the maximal number of solitons is M — 1 for the chain of the length M. 
The isospectral transformation of the spin system depends on the (n — l)-solitonic origin 
and n-solitonic image of the classical states. Finite dimensional similarity operators, called 
the modified Q-operators, of course, do not form a commutative family. Nevertheless these 
isospectrality deformations play the important role in the Functional Bethe Ansatz [|IJ (see 
[|T2| , ||, [13| for the details concerning the method of FBA). Namely, for a special set of the 
solitonic amplitudes the similarity operator becomes a projector into an eigenstate of non- 
diagonal element of the quantum Toda chain monodromy matrix. It was shown in Jl| , but 
we did not give explicit formulas for matrix elements of all finite dimensional operators 
(quantum intertwiner, similarity operators as well as the separating operator). All these is 
the subject of the current paper. 

This paper is organized as follows. First, we recall the formulation of the quantum 
relativistic Toda chain at root of unity and its classical counterpart. Next, we recall the 
form of the auxiliary L-operators (which are the particular case of Bazhanov-Stroganov's 
L-operator [7J), its classical counterpart, the intertwining relations and the isospectrality 
transformations of the Toda transfer matrix. In the third section the matrix elements of 



x in this paper we use the term "r-function" in the sense of bilinear discrete equation with constant 
coefficients, i.e. as a trigonometric limit of the 0-function on a high genus algebraic curve 



4 S. PAKULIAK AND S. SERGEEV 

quantum i?-matrix and modified Q-operator are given in an appropriate basis. In the fourth 
section we recall the parameterization of r-function in application to the relativistic Toda 
chain. In the fifth section we parameterize arguments of the modified Q-operators for the 
whole set of Backlund transformations and so construct explicitly the similarity operator 
transforming the homogeneous chain into the general (M — l)-parametric inhomogeneous 
chain. In the sixth section, we describe an appropriate limiting procedure for the degenerative 
final state and construct the quantum separating operator. Finally, in conclusion, we discuss 
several further applications of the method described. 

1. Formulation of the model 

In this section we recall briefly the subject of the model called the quantum relativistic 
Toda chain at root of unity |], ||. 

1.1. L-operators. Let the chain is formed by M sites with the periodical boundary condi- 
tions, ra-th site of the Toda chain is described by the following local Lax matrix: 

1 _|_ _ u w u \ 

(i) e m (\) = 

\ w m , I 

Here A is the spectral parameter, k is an extra complex parameter, common for all sites, i.e. 
the modulus. Elements u m and w m form the ultra-local Weyl algebra, 



m m i 



(2) u • w = uj w • u 

and u m , w m for different sites commute. Weyl's parameter uj is the primitive root of unity, 

(3) uj = e 2vi ' N , uj 1 ' 2 = e l * /N . 

iV-th powers of the Weyl elements are centers of the algebra. We will deal with the finite 
dimensional unitary representation of the Weyl algebra, i.e. 

where u m and w m are C-numbers, and 

(5) X™ = 1® 1® ... <g) 3C <g)... , Z m = 1 (g) 1 <g> ... ® 2. ®... 

m — th m — th 

place place 



Convenient representation of x and z on the iV-dimensional vector space \a) = \a mod N) 

is 



(6) 



x | a) = | a) u a , z | a) = |a + 1 ) , (a\ 



6, 



a,/3 



Thus x and z are normalized to the unity (x w = z N = 1) N x N dimensional matrices, and 
the Nth powers of the local Weyl elements are C-numbers 



(7) 



N _ N N _ N 



In general, all u m and w m are different, so we deal with the inhomogeneous chain. 

Variables u^ and w^ form the classical counterpart of the quantum relativistic Toda 
chain, and the classical Lax matrix is 



LM N ) = 



( 



\ 



K 



N 



N„.,N 



1 + ^T7 UZW" , -£ 



< \ 



X N 



mm ' 



A^ 



w. 



N 



°/ 



1.2. Transfer matrix and integrability. Ordered product of the quantum L-operators 
d) 



(9) 

and its trace 

(10) 



def 



t(A) = J £ 1 (X)e 2 (X) ■■■ £ M (X) 



( a(A) b(A) ^ 



v c(A) A{X)j 



M 



t(A) = a(A) + d(A) = J2 X ~ kt *- 



k=0 



are the monodromy and the transfer matrices of the model. 

Transfer matrices, defined with the given k and with the given set of u m , w m , m = 1...M, 



:n) 



t(A) = t(X,K;{u m ,w m } m=1 



form the commutative set: V A,/i 



(12) 



t(A, k ; {u m , w m }) , t(//, k ; {u m , w m }) 
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Eq. flT2|) is provided by the intertwining relation for L-operators ([I]) in the auxiliary two- 
dimensional vector spaces with the help of well known six-vertex trigonometric i?-matrix, 
see |I] for the details. 

In the spectral decomposition of the transfer matrix flTTi| ) the utmost operators are 



M 



(13) 



where 



(14) 



t = 1 , t M = k m ] [ (-uj 1/2 u m w m ) Y 



m=l 



M 



y = n (- w " i/2 v-) > yn=i 



m=l 



Further we will need also the monodromy of the classical Lax matrices (§), 

/ A{X N ) B{X N ) \ 



(15) 



and its trace 



T = LxL 2 ■■■ L 



M 



\C{X N ) D{X N ) ) 



(16) 



N\ 



,N\ 



T(X N )=A(X N ) + D(X") . 



M nontrivial coefficients in the decomposition of T(X N ) are involutive with respect to the 
Poisson brackets, associated with the following symplectic form: 



(17) 



M , , 

sr-^ d u m a w m 



m=l 



U, 



W, 



It is the consequence of (0) and ([12|). 

The utmost coefficient in the decomposition of the classical transfer matrix T(X N ), anal- 
ogous to (|T4"D, generates the simple gauge transformation of L-operators ([^). Thus we may 
fix without lost of generality 



M 



M 



;is) 



\{{-^u m ) = \\i 



-V), 



1 



m=l 



m=l 
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1.3. Quantum separation of variables. Eigenstates of off-diagonal element of the mon- 
odromy matrix (|J) b(A) play the important role in the constructing of the spectrum of the 
quantum transfer matrix. This method is known as the method of the Functional Bethe 
Ansatz, or the method of the quantum separation of variables j |T2"| , ^|, [CJ. It is useful to pa- 
rameterize the spectrum of b(A) by its zeros A^, k = 1...M — 1. To deal with the eigenstates 
of b(A) explicitly is not useful. We should better define the set of vectors ({A*;}^ 1 ^), such 
that 



(19) H^lih}^ 1 ,!) = 0, A.GJA,}^ 1 , 



and 



(20) Y|{A fc }ti,7> = \{\ k }fi,l)"\ 

where the integral of motion Y is given by ([14]) . 

The matrix element between an eigenstate of t& and |{A},7) is a product of (M — 1) 
Baxter's functions Q(Xk), taken for the same A^, k = 1...M — 1, as in (|T9|) , see e.g. 0, [13[|. 
For the details concerning the Baxter equation for the quantum relativistic Toda chain at 
root of unity see 0, and later we will recall it. Note here, this Q is a meromorphic function 
on Baxter's curve. In our parameterization it is given by 

(21) x N = 5 N + ^ + k n . 

Convenient Baxter's form x N + y N = 1 + k 2 x N y N may be obtained by the substitution 
x = 1/5, y = S/X, k 2 = -K N . 

Spectra of a, ..., d are described by A, ..., D, see |J3]. For example, 



^22) Yl b K A ) = B ( A 

ngZjv 



N\ 



Eq. ( p2|) allows one to find the set of A&. We will do it later for the homogeneous model. 

The real reason why we are considering the classical counterpart is that the quantum 
separation of variables is associated with the Backlund transformation of the classical model. 
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2. Intertwining relations 

In this section we recall the auxiliary L-operators, quantum i?-matrix and the Darboux 
relation for the classical counterpart, providing the solution of the isospectrality problem, 
construction of modified Q-operators and Backlund transformation. 

2.1. Auxiliary L-operator. Toda L-operators ([!]) are intertwined in the quantum spaces 
with the Bazhanov-Stroganov L-operators 0. This fact is the origin of the relationship 
between the relativistic Toda chain at the Nth root of unity and the iV-state chiral Potts 
model. In general the intertwiner is rather complicated, but there exists a special limit 
of Bazhanov-Stroganov's L-operator such that the intertwiner simplifies. This limit is the 
analogue of the so-called dimer self-trapping L-operators for usual Toda chain, see |TD|, [| ^j 
for the details. We define the auxiliary L-operator as follows: 



(23) 



/ i 1/2 \ 



-cu 1/2 A^u, 1 w w 



CO 



1/2 



uj 1/2 k,wAu, » 



Here A and A^ are two spectral parameters (actually, up to a gauge transformation, £ depends 
on their ratio). k$ is the module as well as k. u^ and w^ form the same Weyl algebra at 
the root of unity (f|) , 



(24) 



u, 



M0 X^ 



W„ 



W<(, 2^ 



In all these notations the subscript stands as the "number" of additional Weyl algebra. 
Classical counterpart of £ is 



/ 



(25) 



-N <t> „,JV 



■it 



N 



N„.,Ny 



l + ^-^w" , -£! + «'> 



\ 



A^ 



A 



N 



\ 



A 



N W <t> 



w 



N 



J 



L-operators (^) are intertwined in their two dimensional auxiliary vector spaces by the 
same six- vertex trigonometric i?-matrix as ([!]). Also there exists the fundamental quantum 



intertwiner for (23): it is the i?-matrix for the chiral Potts model such that two rapidities 
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are fixed to a special singular value. The method describing in this paper may be applied 
directly to the model, defined by (PSf), this is the subject of the forthcoming paper. 



2.2. Intertwining. We are going to write out some quantum intertwining relation for L- 
operators (jl]) as well as for the whole monodromy matrix (§). So we use notations, applicable 
for the recursion in m. Also in this section we will point out parameters u m , w m , (f|), as the 
arguments of L-operators. 

Proposition 1. There exists unique (up to a constant multiplier) N 2 x N 2 matrix Hm,,<f>{^4>) , 
such that R m ,^, £ m and £ obey the modified intertwining relation 

"<t>\^i ""Vi u <j>,mi W cj>,m) ' *-m(^i U mi W m) **-m,<t>\*<t>) = 

(26) 

= ^■m,(j>{^4>) ^mW U mi W m) ' ^(^) *<j>\ U (f>,m+1) W <j>,m+l) ' 

if and only if their classical counterparts of L m and L obey the following Darboux relation: 

(27) Z« m , u#JL«, <) = L(tC, OL(u» m+lJ < m+1 ) , 
Proof: the direct verification. Eq. fl26|) may be rewritten as 

(28) ^(A;u ,w^) -£ m (A;u m ,wJ = £ m (X; u' m7 w' m ) • I^X; uj, wj) , 
where 



"to ""77i """'m,0-'"77i"'""'m,</> ' " m "-'m "-^m, 0^771 "-Sri, <A ' 



V^'-V U m M m ^Hn,cA X m-' r Sn,<A ' W m ^m *Sn,<A Z m-' r S-, 



and 



(30) u^ — u^ m+1 R m ^x^R m ^ , w^ — -u;^ >m+ i R^z^R^ 



10 

Solution of 
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with respect to the primed operators is unique and given by 



(31) 



u„ 



w„ 



u„ 



w. 



K<j> 



u. 



w_w, — a; 1 / 2 A,u, 1 w 



'm vv 



<t> vv * ' 



(\<p + KU m w m - x 



- ,„V2 



w m u^) A^u, 



K 



u m w m (w m u - w 1/2 Aj 



u 



N 



\ N V N 



TVth powers of all primed operators in (31) give 

(32) 

and 

(33) 



.N 



N 



A 6 + K U m W m + W m U <h,rn 



No.N ' W (j>,m+1 



N N 
II, 111 



N \N , ln N v N ' 
K A + W m U 4>,m 



II 



/N 



> JV 

7/ 7/ J 

„/y u '<j>,m i "^m 



/JV 



N N N ln N 
K U m W m "Jfam 



N N N 

K <j> U <f>,m W 4>,m+1 



Eqs. Q32|) and ( |33"D are the exact and unique solution of (|27|). Therefore eq. (|27|) is the 
consistency condition for (p6| ). Next, eq. ( p6|) is a set of linear equations for the matrix 
elements of R m ,^, this set is defined appropriately, and this provides the uniqueness of R m ,<^. 
The local transformation 



u m, w m, u <j>,m, w (j>,m h ^ y u m i W mi U <j>,m+1, w <p,rn+l i 



(34) 

given by eqs. (|3~2"D and (|33|), is called usually as the Darboux transformation for the classical 
relativistic Toda chain. Eqs. ( |3"2D and ( |33"D define the mapping ( |34"D up to iV-th roots of 
unity. These roots are the additional discrete parameters of the transformation fl3"4|). Note, 
matrix ~R, m ,<t> is unique if all these roots are fixed. 

2.3. Q-transformation. Relations ( p6|) and ( p7]) may be iterated for the whole chain. The 
functional counterpart gives 

(35) £(«&,<!)?(« «£}£*) = T({u^,w^}^ =1 )L(u^ M+1 ,w^ M+1 ) , 



where the mapping u^^.w^^ i— ► tt^Af+i, iu^m+i is M-th iteration of (|34|). 
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For the periodic chain the cyclic boundary conditions for the recursion fl32| ) are to be 
imposed, 

(36) U^m+1 = U^i , W^m+1 = Ul^i . 



Now suppose ( 52|) and ([5B) are solved, i.e. M^ m , W0 jm are parameterized in the terms of 
u m , w m , m = 1...M, and some extra parameters, possible degrees of freedom of CJ32|) and (|36|) 
(e.g. A^). Then ([^) defines in general the transformation Q^ 



(37) Q<f, : {-u m ,w m } m=1 i— >• {u m ,w m \ m=1 

The transfer matrices of two sets, {u m , w m } and {u' m , w' m }, have the same spectrum, because 
of there exists N M x N M matrix Q^, nondegenerative in general, 

(38) Q^ = tr<£ (Ri,^ R2,</) • • • Ra/,</>) , 
such that 

(39) t{X;{u m ,w m }) Q^ = Q^ t(A ; {u' m , «4}) . 

Subscript of Q^, stands as the reminder for the parameters, arising in the solution of 
recursion, including at least the spectral parameter A^. 

Consider now the repeated application of the transformations Q^, eq. (|37|), 

(40) {u m ,wj % {u' m ,w' m } ?t 2 ... & {«£>,«;<?} , 
such that the set of isospectral quantum transfer matrices 

(41) tW(A) = t(A,«;{ri»),u)W}li) 
has arisen. Sequence fl40|) defines the transformation /C, 

(42) /C^ : {u m ee u£\ w m ee ui 0) }™=i ^ {«£>, ^ n) }^=i , 
with the finite dimensional counterpart 

(43) K<»> = QgQg.-.Qj?, 
where 

(44) Q^ = tr^ n (R^ R 2 ^ n . . . Rm,^ j 
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makes 

(45) t^WQj? = Q^t^A), 
and 

(46) t^(A) KW = K("» t^(A) . 

3. Matrix R 

In this section we construct explicitly the finite dimensional matrix R mi< £, obeying fl2~6|) , 
in the basis (|6]). But first we have to introduce several notations. 

3.1. ^-function. Let p be a point on the Fermat curve JF 

(47) p = f (x,y) eT <£> x N + y N = 1 . 

Very useful function on the Fermat curve is w p (n), p G J-, n G Z^r, defined as follows: 

(48) W ; {n) , = y , Wp (0) = 1 . 

w p (n — 1) 1 — xuj n 

Function w p (n) has a lot of remarkable properties, see the appendix of ref. |15| for an 



introduction into cu-hypergeometry. In this paper it is necessary to mention just a couple of 
properties of w-function. Let O be the following automorphism of the Fermat curve: 

(49) p = (x,y) O Op = {uj^x' 1 ^uj^^x^y) . 
Then 

( 50 ) w 0p {n) = —— — , 

®(n)Wp(— n) 

where 

(51) $(n) = (~) n u n2 / 2 . 

In the subsequent sections we will use implicitly two simple automorphisms else: 

y 

(52) W( x>uy )(n) = u n w^ y ){n) , W( XtV ){n + 1) = — W{ux, v )\n) ■ 

Define also three special points on the Fermat curve: 

(53) g = (0, 1) , goo = Oq , gi = (a; -1 , 0) . 
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Then 



(54) 



w qo {n) = 1 , w qco (n) 



1 



1 



$(n) ' w qi (n) 



Jnfi 



3.2. Matrix elements of R m ^. Consider the N 2 x N 2 matrix 
Rm,4>(PiiP2->P3) with the following matrix elements: 



(a m , a^lUmrflPm, (3$) 



(55) 
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{cLm-MPt Wpi((*<l> ~ <*m) WptJP,!, - /3 m ) 

w P3 (fy - a m ) 
Here Px,P2,P3 are three points on the Fermat curve, such that 



'atStPt 



<£,pm 



(56) 



Xi x 2 = x 3 . 



Eq. fl5"6|) and the spin structure of (p5|) provide the dependence of fl5"5"D on two continuous 
parameters, say xi and x^, and on two discrete parameters, say the phase of y\ and the 
phase of y 3 . 



Proposition 2. Matrix R*m,$(pi,P2,P3), whose matrix elements (51.) are given in the basis 
(|6[), makes the following mapping: 



(57) 



-i 



"^m^^-m^m 



-TV-™ aZ_,JX„ j, — Z„,Zj, UJ 



m,q> m m,4 



"-m,0 X *V,<^ 



*Sn,<£ z </, ■"'m,^ 



</2 



ggj/l -1 
X W 2 ^ '< 



_i a;i2/ 3 _! 



u;x 3 x m -w 



2/1 



2/3 



-1 _ 2*J/l -l -l 

X X m. W Z m X m 

2/2 V X\V2 



Proof: direct verification with the help of (|52f) . 
Compare ( |57l) with ([Slf ). Obviously, R m ,0 solves (|26|) if 



2/3 r 

2/i ' 



(58) 



Si. 



CJ 



-1/2 ^</>,m ^3,m 2/l,m 

K U m X\ rn y 3 ^ m 



U) 



-1/2 *V 



"'(f/,m"^m 
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and 
(59) 
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u, 



K 



U, 



w„ 



w„,w. 



J/2,i 



m^0,m 



U 4>,m+1 ' = ^ x 3,m U m , W0 iTO +l 



2/3,, 



K M m 2/3, m 



^</i u 4>, m 1j2,m 

Eqs. (|58|) and ([59]) allow one to parameterize R^, in ([44]) in terms of u m) w m) .... So for 
any appropriate set of {u m ,Wm }, m = 1...M and /c = 0...n, see formula (^0|), we may 
construct the corresponding K> n ) explicitly via the parameterization (|5"8| , |59"D , form of the 
matrix elements (|55|), and the formulas ( [44"| , [4"3|) . 

Further we will consider the sequences like (f40|) for the homogeneous initial state 
{u m = u ,w m = u> }m=n an d parameterize {um \wm } in the terms of solitonic Hirota- 
type expressions. To do it, several notations are to be introduced. 



4. Rational 6-function 

In this section we introduce several useful functions and notations. The main purpose is to 
introduce the form of solitonic r-function for the classical relativistic Toda chain. The reader 
may find the details concerning the corresponding classical integrable model, reduction from 
2DToda lattice hierarchy and so on in Ref. [{H|| . 

Fix any sequence of (pk from the set 3m : 

def ( 7T 2n 



(60) 

Introduce the functions 



G S 



M 



{ 



M' M 



(M-I)tt 
M 



} 



(61) 



cos z 



+ K N + COS( 



< a; = e - ^ L/co&cj) + k n + cos^) , 



a^ = a^ a; . 

Expressions for A and A* parameterize the curve 
(62) AA* 



A 



A + A* + k n in the terms of e 2i0 = — 

A* 
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Rational B-functions are to be denned recursively as 

e (0) = 1 

(63) e£> = 1 - h e 2im<t)1 , 

eg) = 1 - f\e 2im ^ - f 2 e 2im ^ + rfi,2/i/ 2 e 2im(01+02) , 



and so on, 

(64) ew = e^\{f k , 4> k y k z\) - /ne^e^^d^, <j> k } n k z\) , 

where all 0^ are different and belong to the set (^Dj). The phase shift dj tk = d^.^ k is given 
by (see flBJ) 

(65) d - (A * ~ V) ( A ~ A l') 

Arguments of 6™ are the set of 4> k G #m an d corresponding amplitudes f k , k = l...n. 
Maximal Q-i unction is 

(66) e m d = f e^\{f k }iLS) . 

^A/-set, (|S0|), provides the M-periodicity of B-functions, Q m+M = m . 
Introduce next a couple of useful functions 

A* A , - f 

(67) s&) = / ^J— i , 

and 

A* Al - A . 

(68) a*, = s,(a;,) = -± -I -J . 

Necessary is to consider the re-parameterization of the amplitudes f k in the special form: 
for any n, k, < n < M, 1 < k < n let 



(69) ft\{g}) = g k f[ 



°fe,7 



J ' 






where s^j is given by (0). Denote the corresponding 0-functions as follows: 

(70) &£>M) = eW({/ fe w ({^}),0 fe }Li) 
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In advance, we will always imply the free amplitudes g^, k = 1, ..., as the arguments of 
6-function, (^,[TtJ). 

Particular case when all gu = Sfc(£) is also important: let 

(7i) eS?(£) = eW({y* = Sjfc (e)}Li) 

This particular case has the following property: 
Proposition 3. For any n 

(72) e«( 9 = n|^ n #^#. 

fc=l k S l<j<k<n 3 k 

and for m = 0...n 

(T3) eg)(0 = (-^) (n " m )(n ~ m - 1)/2 (^^J ©i n) (£) 

//£ = 1, eg. ( f75^ zs wa/zd a/so /or m = — 1. 

5. Parameterization of K 

5.1. Sequence of the Backlund transformations. Turn now to the original quantum 
chain. Let the initial parameters u m ,w m are homogeneous: 

(74) u m = -u~ 1/2 , w m = -1 . 



It corresponds to (|i8|). 

For the sequence of (pk € $m fixed, and for generic set of corresponding o^, k = 1...M — 1, 
let (see eqs. (p9| ) and (f70| ) and the note right after flTDD) 



.AT „^,, ,„ , ,„^\AT 



(75) ffl = e£>({**}Li) , (4 n) ) = ©L n) ({^(i)}Li) • 

The phases of Tm and #m are arbitrary. 

Proposition 4. Consider the set of (M — 1) transformations fiSlj ), such that the initial state 



is the homogeneous one (\74), and the mapping Q^ n is derived for 
(76) A£ = A. , 
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where sequence of (fi n G 5a/ is fixed, and function A^ is given by the last formula of ffil\). Let 
besides 



(77) 



.N 



Then the sequence fc{(\ ) of {um\wm } may 6e parameterized as follows: 






(78) ut ] = -u~ l/2 ^-, w™ = --£-, ; mez 



(ft) ' m /)(") 



^Ai 



'm ^m-l 

where Tm , Om are given by ff7b\ ) with arbitrary gi-.gM-i- 

See for the sketch of the proof of this proposition. Eq. ( fFTh needs a comment. Arbitrary 
value of K0 n corresponds to a gauge transformation of 

(.'"J rm > W m lm=l l—> ' l CM m > C M m lm=l 

with some c proportional to «,*„, while the structure of r and 9 is not changed. So we impose 
Q77| ) and obtain ( [78] ) without lost of generality. 

5.2. Parameterization of modified Q-operators. Let now for the shortness A^ n = A n 
are given by (|7|), and let us fix further the roots for S n and c n , 1 < n < M : 

(80) *n=V> «*. = f > C ft=%^" 

Auxiliary values in (|^,^) are given by 

(ft) 

/o-i \ -1/2 r T m-1 -1/2 

(81) «^, m = -w ' n ^^ , W^ tt , m+1 = W 7 

.(ft) 

x -lllclljl 1UCD, CllljCllllg, IjU \tV-i 

W) R (n) ■ n {n) n {n) n {n) 

3Z 7 m,0 n • x v l,m; i^.m) F3,mi 



r ( n ) JC /d(«-i)_(«) 

' m— 1 —1/2 ft m ' rn 

(7^ , «W+i = a; — 

<rn " " 'm "m 

R™, -matrices, entering to (fB|),have the arguments 



given by 



(83) 



x 



Un ' m '■>), — '! (rj.l _1 r- ' -r^ ' m 



(ft) -1/2 £ft 'm 'm-1 (n) -If 

l,m ~ W K (n-1) (ft) ' X 3,m " ^ °« (n-1) (n) 

'm—1 m m— 1 'm+1 



„(») JC fl(ft-l)-(ft) 7 / n ) X fl («-l)_(n) 

^3,m _ ,-1/2 "ft ^m 'm-1 #l,m _ 1/2 ^ft "m-1 'm+1 

(n) ' A ' \ n (n— l)/i(n) ' (n) ^' ^.X /i(ft— 1) (ft) 

t>2,m m—1 m t>S,m m m 
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and with fl56|) implied, 

m—1 ra+1 

The trace of quantum monodromy (|44] ) may be calculated explicitly in the basis (§), and 
the answer is 

nWM {(3 m - a m ) w („) (/3 m+ i - /3 m ) 
^HW/Wl _Pm. ?v, 

WW /W - OC m ) 

m£Z M P3,m 

Operator K, (f£3"[), calculated as the product of all (M — 1) Q^-opeiators, is thus defined 
explicitly. 

Proposition 5. Operator K({g}) ; gwen &?/ (|J|j ; does not depend on the ordering of(gk, <pk) 

(up to a multiplier, arising in general when one takes arbitrary phases for Tm and 9m ). 

Tin 
Thus without lost of generality one may regard n = — . 

Thus for the generic set of amplitudes gt constructed is the similarity operator K({o}), 
transforming the homogeneous initial transfer matrix into most generic isospectral one. 

5.3. Arbitrary value of A^. Previously we considered the set of 0& G 5m- Actually, all the 
calculations in the classical relativistic Toda chain are based on (p^). Therefore, in general, 
we may consider a generic sequence of complex numbers 0^, k G Z + . For any generic complex 
one has to restore A^ and 8$ via (|6l|), and take into account (J30|) . It corresponds to the 
Q operator for the generic value of A^,. Then all the formulas for the parameterization of 
transfer matrix in the terms of r-functions (|78|j75| , [70|) as well as the parameterization of Q- 
operator (|33| , |85D are valid if the amplitudes o^, corresponding to those G" $m an d — G" $m 



are zeros, and in the case when G $m or —0 G $m the amplitudes are arbitrary. Important 
is that in the calculation of the amplitudes ( |69D all 0s are to be taken into account. Also, 
if in the generic sequence of 0^ it happens 0^ = m or fc = — m , then 6-function is to be 
understood as the corresponding residue of formal singular expression (|70|) . 

For example, usual Q-operators for the homogeneous initial state corresponds to 0i = 
^ -3a/) 9\ — 0. The cases when 0^ = ±0 m correspond to the annihilation of the solitons. 

Nevertheless, considering operator K, (H3), we will always imply 0^ G $m ■ 
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5.4. Baxter equation. Consider initial homogeneous transfer matrix t(A) and correspond- 
ing Q operators with arbitrary value of A^. As it was mentioned in the previous subsection, 
corresponding functional counterpart of Q-operator is trivial, i.e. such Q-operators commute 
with the transfer matrices and form the commutative family. Such Q-operators are usual 
Baxter's Q-operators, obeying Baxter's equation. 

Explicit form of simple Q-operator is given by ( [35"1) with the homogeneous parameteriza- 
tion 

(86) p 1>m = pi , p 2 , m = P2 , P3,m = P3 

with 



-l/2"0 -1/2 -1 

X\ = CO ' — - , X 2 = LO ' K , X 3 = CO 

K 

(87) 

- l 'Hs w, ^ 1/2 A, 



2/3 _ u ' 0<t> 2/i to 
i 



V2 A^ y 3 kS^ 
where 6$, X^ and c^ are given by flSOlfTB"! ) and (|61|), but with generic value of 0. Using 



we may regard 

(88) Q = Q(A^). 

Several simple properties of simple Q-operator may be derived with the help of the matrix 



elements of (JS5|) and fl87|) , and with the help of (]52|) : 



XQ(A^,^)X" 1 = Q(uj- 1 X^5 4 



(89) ZQ(A^,^)X 




YQ(A ,^) = 


f \ 1 fi \ M 


where 






M M 


(90) 


X = ] [ x m , Z=][z m , 




m=l m=l 


and Y is given by (fbl|). 





20 S. PAKULIAK AND S. SERGEEV 

Let us derive now Baxter's equation for general inhomogeneous chain. To do it, one has 
to consider the degeneration point of £, (B3|): 

1 - U^K^Wj, 

Therefore eq. (|26|) in the degeneration point may be rewritten in the following forms: 



(91) 



^>(A</,, -"V) 



-CO 



1/2 AtV 



(92) 



and 



(93) 



1 ' ~ U \ X < 



I ' <- m (A^,M m ,w m J • rt r 



R 



(i) 

m,<j) 



1 ,.l/2 U $,m+l 



/ l/2 U <P,m+l \ 

U A, ** 



£m{^<t>\U m ,w' m ) ■ R r 



( 0,1/2^ \ 



R 



(2) 



where 



(94) 



R 



(i) 

m,<j) 



U, 



u 4>,m+l 



■^■■m "rn-rh X 



-1 



R 



(2) 



1 12 ^tb.m+l W m 
W x Z m *W X^ 



Let further Q \ and Q\, are the traces of the monodromies of R^ ', and R^ ' , (Zm boundary- 
condition for Ufcm is taken into account again - the subscript <fi reminds this). Then from 
(|52]) and ( |9"3"D it follows the Baxter equation in its operator form: 

(95) t(A ) Q^ = Q t'(A^) = q£> + Q< 2) . 



Using the first relation of (57), one may obtain 

(96) qW = JJJ^ X QX- 1 , Q(2) = -Q^i/2!f^!^ 



"</>,'■ 



A„ 



ZQX 



where X and Z are given by (|90|) 



21 

Turn now to the homogeneous chain. Using (|89|), we obtain 

/ r \ m 

(97) t(A*)Q(A ) = <^ M Q(uT%) + I -^^f ) Y QM^ > 

where (^-argument of all Q-s remains unchanged, and t and Q may be diagonalized simul- 
taneously. Let t(X) and <ft(A) be the eigenvalues of t and Q for the same eigenvector, then 
(|9"7| ) provides the functional equation 

/ e\ m 

(98) t(X)q t (X) = 5- M qt{^ 1 ^) + f -^ 1/2 j J ^q t {uX) , 

where uo 1 is the eigenvalue of Y, and qt is the meromorphic function on the curve (J2ip . 



Detailed investigations of (pq) show that for generic A any solution of (98) such that £(A) is 
a polynomial and (ft (A, <5) is a meromorphic function on the curve fl21|), gives the eigenvalue 
oft(A) andQ(A). 

6. Separating operator 

6.1. Spectrum of b. For the homogeneous chain ([74]) the value of B(X N ), fll5|), may be 
calculated immediately: 

M-l 



(99) B(A") = (") 7V - 1 ^ II f 1 



A, 



\N 11 V X N 

fc=l v 

where A^ is given by the last equation of ( |6Tj ) and 0^ is exactly the set $m, (|60|). This defines 
the parameterization of the vectors ({Afc}^ 1 ^), see (|HJ|). 
In [|T| it was established the following fact: 

Proposition 6. Let A^ = A^ fc are defined by ffTBJ and the phases of A& are fixed. Operator 
K{ A }, corresponding to the case when all g^ = 1 , k = 1...M — 1 in eqs. ( [7o[ , |77| j, makes the 
quantum separation of variables: 

(100) b(A fc )K {A} = 0, k = l...M-l. 

This means 

(ioi) k {a} = J2 \{^}k=fn)(x,\, 



22 S. PAKULIAK AND S. SERGEEV 

where x-y are some vectors, such that (x 7 | Y = u; 7 (x 7 |, where Y is given by (|14|). Parame- 
terization 

(102) {r^) N = eW(o) , (*£>)* = eW(i) 

is not appropriate for (|S3D, because of ri = for m = 0...n — 1, and Um are ambiguous. 



Thus one needs to define in addition a limiting procedure for ( |102| ). In general, the ambiguity 
in (|78| ) corresponds to the arbitrariness of the right eigenstates x-y °f K{a}. 

6.2. Limiting procedure. In this section we describe an appropriate limiting procedure. 
Most simple K{ A } appears when we choose u m = for m = 1...M — 1. Namely, consider 
the set of infinitely small numbers, 

(103) e k i-> , k = 1...M-1 , 

such that any ratio e n /e m ^ 1 is finite, and 



j{M-l) = , .-1/2 £ M 

so that 



(104) u^-V = -oj-V 2 - &AJ : m , m = l...M-l, 

Ol....() A f-l 



(105) 4f- 1} = ,^-1/2 (^i-^-i)^ 1 

£l...£M-l 



Recursion ( |32"| , |33"D with 

(106) «*, = f 

On 

implies 



(107) «&) = -u- 1 ' 2 ^=£1 + o( £l ) , m = l...n 

Oi...0 n 

while for m — n+ 1...M — 1 all Um are regular. This means 

win as explicit functions of Sk are rather complicated, and speaking truly, we may say 
nothing about t( M ~ x ) in the meantime. But in the limit e^ = the following formulas may 
be chosen: 

(n) n 

(wo) * = n c t , 

° n fc=i 



23 

where c^ are given by fl50|) and 



(110) 9$ = ( W V2 K )(n-m)(«-m-l)/2. 






(5 1 ...5 n ) n - m ' 
On iV-th powers these formulas follow from proposition (|3|). Eq. ( |110| ) implies in part 

w {n) 
(111) -Zp± = co 1/2 K, m = 0...n-l. 

( n ) 

Win 

Now in the limit Sk i— > the following parameterization of the ambiguous and singular 
points arises (notation for the singular points are described by fl53| , |54D ): 

mm — l...n — 1: 

(112) 



where x^m = 
m = n: 


(n) 
P3,m ~ 


= ?1 , 


(n) 


= o(p?i) 


, il/2if r 

If rv tn— m 





(H3) pi5 = go , p£! = g^ Vn 

and p 3 "„ are regular, except 

(H4) pi M M '-l = go 



loo 



In) 

m = n + 1...M — 2 : All p) ^ in this region are regular. 



],m 



m = M — 1 : Pil[_i are regular except 



(us) p[ M m-1 = go , 



and 

(116) P^m-i = Ptli-i = g°o V n . 

• m = M : 

/ 1 , 7 \ (n) (n) (n) w 

(11 7 ) Pi,Af = g°o , p2,m = Psif = go Vn. 
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Substituting these expressions into (|85|), one obtains the following form of n-th modified 
Q-matrix, n^M-1: 

n-l 



m=l 



^m + l ' 



1181 



$(«„) tfppW(«n - Pn+l) X 






$(«M- 



-T W (n) (0M-1 - OiM-l) 



and for the last (M — l)th one 

1 Af-2 

(119) (ciq^- 1 ^) = — — - — -— n ^.*w • 

Here $ is given by (|5lD . Explicit form of the modified Q-operators ( 118 ) and (|119|) allows 
one to prove the following 

Proposition 7. 

(120) K {A} (z m - z m+1 ) = , m e Z A/ . 
So, in our limiting procedure x is simple: 

(121) (x 7 |«) = X 7 («) = X 7 (« mod #)> 
where 

(122) at/ £ 

and 

(123) X 7 (5 + M) = ( _ w - 1/rw7 - s 

X 7 («) 



a, 
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7. Discussion 

In this paper we investigated the simplest integrable model, associated with the local 
Weyl algebra at the root of unity. All such models always contain the classical discrete 
dynamic of parameters. Nontrivial solution of the classical counterpart provides the solution 
of the isospectrality problem of the finite dimensional counterpart. Well known result by 
Sklyanin, Kuznetsov at al was that in the classical limit of the usual Toda chain (and many 
other models) Baxter's quantum Q-operator corresponds to the Backlund transformation 



of the classical system, see e.g. M, [LOl 11, 12, [Lq] etc. In our case we have the Backlund 



transformation of the classical counterpart and modified Q-operator in the quantum space 
simultaneously. Unusual is that solving the quantum isospectrality problem, we miss the 
commutativity of the modified Q-operators. 

Nevertheless, our main result, we think, is the explicit construction of (M — l)-parametric 
family of quantum inhomogeneous transfer matrices with the same spectrum as the initial 
homogeneous one, and the explicit construction of the corresponding similarity operator 
(|43|). We hope, the solution of the isospectrality problem will help to solve the model with 
arbitrary N explicitly. 

As one particular application of it we have obtained the quantum separation of variables. 
Previously there was a hypothesis, formulated for the usual quantum Toda chain, that the 
product of operators Q, taken in the special points, is related to the quantum Functional 
Bethe Ansatz. In this paper we have established it explicitly, but for the product of modified 
operators Q. Below some explanations are given. 

Let the solitonic amplitudes gt are in general position. Let \W t ) be the complete set of 
eigenvectors of t^ n ^(A), and (\&£ | be the corresponding set of co- vectors, 

(124) &\X) = J2 \*l n) )t{\)(9l n) \. 

t 

Then 

(125) QjJ = J2 \^ n - 1] ) Qt(Xn) (^ n) \ , 

t 

where (ft (A) is evidently the eigenvalues of simple Q-operators, as in (0). Note, the functional 
Baxter equation exists also and for arbitrary inhomogeneous chain, but to obtain it from 



95|) one has to use the decomposition like ([12 
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With ( |125| ) taken into account, K-operator for general set of g k is 

(126) K = E l*f >fc(AO...ft(AM-i) (*r _1) l , 

t 

In the limit when all ^ = 1, the sets (W* ? \ become degenerative, e.g. 

(127) (*i n) |(zi-z M ) = Vt. 

For us it is still mysterious what happens with the set |\&t ) and co-set (f{ when g^ — 1 
- this is the subject of separate investigation. But nevertheless, in [[jj and in this paper it 
is proven that K is given by ( |1 26| ) and simultaneously by (|101|) . Thus Sklyanin's formula 
appears explicitly: 

Af-l 

(128) (^ 0) \{h}t-i\l) = const J] q t (X k ) . 

fc=i 

Well known is that to solve Baxter's equation for iV > 2 is quite hopeless (pure fermionic case 
N = 2 is trivial, flffif) may be solved in one string). Only way is to avoid Baxter's equation 



for high genus curve. So, we think, the attention is to be paid to the special cases of the 
isospectral transfer matrices like g k = 1, when several simplifications for the eigenstates are 
expected. 

Note in conclusion, this method may be applied to any model, based on the local Weyl 
algebra. Mentioned are to be the chiral Potts model || [7| and the Zamolodchikov-Bazhanov- 



Baxter model in the vertex formulation [15|. To be exact, all two dimensional integrable 
models with the local Weyl algebra are particular cases of the general three dimensional 
model, and their classical counterparts are known ||. 
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